Localization and One-Parameter Scaling in Hydrogenated Graphene 
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We report a metal-insulator transition in disordered graphene with low coverages of hydrogen 
atoms. Hydrogen interacting with graphene creates short-range disorder and localizes states near 
the neutrality point. The energy range of localization grows with increasing of H concentration. Cal- 
culations show that the conductances through low-energy propagating channels decay exponentially 
with sample size and are well fitted by one-parameter scaling function, similar to a disorder-driven 
metal-insulator transition in 2-dimensional disordered systems. 
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Graphene, a single layer of graphite, has its unique 
electronic structure with a zero gap and quasiparti- 
cles described by massless Dirac fermionsjii^ The lin- 
ear dispersion relation near two inequivalent Brillouin 
zone corners leads to very unusual transport phenom- 
ena. Dirac fermions can perfectly transmit through a 
potential barrier at normal incidence due to the ab- 
sence of backscatteringj 3 . so called the Klein paradox. 4 
The metallic conduction of graphene is robust against 
long-range disorder, whose potential varies slowly on the 
scale of the interatomic distance^! 6 - Thus, Dirac fermions 
cannot be generally localized by long-range disorder, in 
contrast to the Anderson localization theory^^ while 
a metal-insulator transition of Kosterlitz-Thouless type 
was recently reported in graphene with strong long-range 
impurities^ On the other hand, random short-range dis- 
order induces the intervalley scattering between the two 
valleys, eventually leading to localization ! 10 i 11 

The chemical reaction of atoms and molecules with 
graphene can give rise to a short-range disorder poten- 
tial. Several experiments reported that graphene under- 
goes a metal-insulator transition by dosing with atomic 
hydroge n 12 ' 13 or molecules such as NO^J^ At high cov- 
erages of adsorbates, the electronic structure of graphene 
may be significantly modified due to the change in hy- 
bridization from sp 2 to sp 3 . According to first-principles 
calculations^ when all C atoms react with hydrogen, 
graphene turns into a new insulating material with the 
band gap of about 3.5 eV, known as graphane. Other 
theoretical studies showed that chemisorbed molecules 
such as H and OH suppress the conductivity on one side 
of the Dirac point and derive the system further towards 
the localized state at higher adsorbate concentrations. 16 
As the insulating behavior of hydrogenated graphene was 
observed at substantially low doses,— there is other pos- 
sibility that short-range disorder by hydrogen induces a 
localized insulating stated 

In this paper we perform numerical calculations to in- 
vestigate the localization behavior of disordered graphene 
by hydrogenation using a simple tight-binding model. To 



exclude the formation of a band insulating graphene, we 
consider low coverages of atomic hydrogen. We find that 
conductances in a narrow range of energies near the Dirac 
point are well described by one-parameter scaling func- 
tion, exhibiting a metal-insulator transition. 

Our calculations are performed using a combined ap- 
proach of the density functional theory and the tight- 
binding (TB) method. As the electron conduction takes 
place by hopping along the C ir orbitals, we consider a 
single-band TB Hamiltonian to describe interactions be- 
tween graphene and hydrogen, 
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Here 7 (= 2.6 eV) is the hopping integral between the 
nearest-neighbor C it orbitals, Ci (Cj) is the annihila- 
tion (creation) operator on the Ith site of graphene lat- 
tice, is the coupling strength between the C and H 
orbitals, e# is the H on-site energy, d n (d^ ) is the anni- 
hilation (creation) operator on the adsorbate site, and 
p n is the host site bonded to the H atom. The pa- 
rameters, 7ff (= 5.72 eV) and e# (= eV), are de- 
termined by fitting to the first-principles band struc- 
ture of hydrogenated graphene. In first-principles cal- 
culations, we use the generalized gradient approxima- 
tion (GGA)ii for the exchange-correlation potential and 
ultrasoft pseudopotentials 18 for the ionic potentials, as 
implanted in the VASP codeJ^ The wave functions are 
expanded in plane waves with an energy cutoff of 400 
eV. We test various hexagonal supercells with different 
sizes, which contain up to four H atoms, and find good 
agreements between the TB and GGA band structures, 
without including the second nearest-neighbor hopping 
parameter in the TB model. 

To describe the transport properties of hydrogenated 
graphene, we set up a device model such that disor- 
dered graphene is sandwiched between two semi-infinite 
graphene electrodes, as shown in Fig. 1. In the de- 
vice region, a square shaped sample with the length L 
is assumed, with periodic boundary conditions imposed 
in the transverse direction to remove the effect of the 
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FIG. 1: (Color online) A device model for disordered graphene 
with random hydrogen adsorbates between two semi-infinite 
graphene electrodes. 



edge states in sample. For various coverages, hydrogen 
atoms are placed on one side of graphene because the TB 
approach cannot distinguish the difference between both 
sides. For large systems, we use the recursive Green's 
function method^ to calculate the density of states and 
the two-terminal conductance, <?l = 2Tr(tt^), where t is 
the transmission matrix and the factor 2 accounts for 
spin degeneracy^ 

A single H atom on graphene breaks the sublattice 
symmetry, opening a band gap and creating a localized 
level at the Dirac point. The wave function amplitude of 
the localized state is zero in the same sublattice as the C 
atom bonded to hydrogen, whereas it decays rapidly with 
distance from the adsorbate in the opposite sublattice, 
similar to an ideal C-vacancy^ This similar behavior re- 
sults from the formation of a covalent C-H bond involving 
the p z orbital of the C site where H is adsorbed, effec- 
tively removing one C atom from the lattice. If a small 
amount of H is adsorbed on graphene, the formation of 
localized levels strongly depends on the type of adsorp- 
tion sites due to their unique wave function amplitudes. 
The densities of states and two-terminal conductances 
are compared for graphene samples with various H den- 
sities (riH up to 10 %) in Fig. 2, where «h is defined as 
the ratio of the number of adsorbates to the total number 
of the C atoms in graphene. We consider two different 
configurations for the positions of adsorbates, which are 
in the same sublattice (A or B) or in both the A and 
B sublattices with equal amounts. For each concentra- 
tion, the adsorbate sites are chosen at random, with the 
H atoms on top of the host atoms. In the former con- 
figuration, all the localized states induced by adsorbates 
are positioned at the same energy, E = [Fig. 2(a)], 
because these states are decoupled in the same sublat- 
tice. The van Hove singularities which appear at E = ±7 
in clean graphene become softened as nn increases. Al- 
though the selective dilution of adsorbates is unlikely to 
occur, it is interesting to note that adsorbates develop the 
gap opening. The energy gap has a tendency to increase 
with increasing tt-h-— In the gap region, the densities of 
states are zero, except for strong peaks by the localized 
states, and conductances are severely suppressed [Fig. 
2(c)], indicating that hydrogenated graphene becomes a 
band insulator. 
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FIG. 2: (Color online) The densities of states for graphene 
samples with hydrogen atoms randomly distributed on (a) the 
same sublattice sites and (b) on both the A and B sublattice 
sites. The normalized conductances, <?L/<?E lax > where <?™ ax is 
the maximum conductance of clean graphene at E = ±7, are 
plotted as a function of energy in (c) and (d) for samples 
(a) and (b), respectively. The insets show the logarithmic 
plots of gi^/g^ 3 *- Lines with open circles, triangles, squares, 
and filled circles correspond to the H concentrations of 0.5, 
1, 5, and 10 %, respectively, whereas plain black lines are for 
pristine graphene. 



When adsorbates are random in both the A and B sub- 
lattices, the localized states have nonzero wave function 
amplitudes at the opposite sublattice sites. Due to the 
level splitting by interactions, the gap opening is sup- 
pressed, with several sharp peaks superimposed on the 
finite density of states near the Dirac point. This feature 
is very similar to the case of random vacancy defects^ 
In the region of high energies, the variation of normal- 
ized conductances with energy is similar to that for the 
H adsorption in the same sublattice, regardless of the H 
concentration. The logarithmic plot of normalized con- 
ductances shows clear differences, especially in the low 
energy region, where the finite density of states is formed. 
The logarithmic conductances severely fluctuate, indicat- 
ing a signature of localization, as shown in Fig 2(d). The 
energy range of fluctuating conductances grows with in- 
creasing of the H concentration. In addition, although 
conductances are larger than those for adsorbates in only 
one sublattice, they are much suppressed due to scatter- 
ing between the localized states. 

To see more precisely the localization behavior of 
low energy states, we examine the hypothesis of one- 
parameter scaling which has been widely used in 2- 
dimensional (2D) disordered electronic systems.— In the 
scaling theory of localization, we consider the intrinsic 
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FIG. 3: (Color online) The averaged conductances are plotted 
as a function of sample size (L) for the H concentrations of 
(a) 0.5, (b) 1, (c) 5, and (d) 10 %. In (c) and (d), the insets 
show the localization lengths (£) for different energies. 



conductance g^ which is given by the relation, i = 
— — 7rir, where N c is the number of channels at the 
energy E and j^j- is the contact resistance. Using the 
dimensionless conductance g, the scaling function (/3) is 
defined asZ&& 



d{\ng) 
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Here (...) denotes the ensemble average of lng over config- 
urations chosen for the random distribution of adsorbates 
in both the A and B sublttices of graphene sample with 
the size L. We test various system sizes up to L — 44 
nm, which contains 7.5 xlO 4 sites. For each adsorbate 
concentration in a given sample size, we use 300 to 1000 
configurations, which ensure the numerical convergence 
in the average. 

Figure 3 shows the variation of (lng) with the sample 
size for different energies to within 0.8 eV from the neu- 
trality point, in which the density of states is greatly af- 
fected by adsorbates. When adsorbate densities are low, 
fin < 0.5 %, conductances are slightly fluctuated due to 
small sample sizes for L < 20 nm [Figs. 3(a) and 3(b)]. 
For Tin = 0.5 %, (lng) remains nearly constant with in- 
creasing L, regardless of the energy E in the vicinity of 
the Dirac point. This result is consistent with the fact 
that the density of states is not significantly affected at 
low densities of adsorbates, with only sharp peaks at the 
Dirac point [Fig. 2(b)]. As hh increases to 1.0 %, (lng) 
with E = —0.2 eV starts to decrease with L, while no 
change occurs for higher energy channels. The decreasing 
behavior of (lng) with E = —0.2 eV becomes significant 
for high adsorbate densities of 5 and 10 %, as shown in 
Figs. 3(c) and 3(d). It is interesting to note that the 
energy range, in which (lag) decreases with L, increases 
with increasing Tin- For njj = 10 %, we find the lin- 
early decreasing behavior of (lng) for all the propagating 
channels down to E = —0.8 eV, with different slopes. As 



the conductance decays exponentially with respect to the 
sample size, the slope represents the inverse of localiza- 
tion length (£). The localization lengths are estimated to 
be in the range of 2 to 7 nm for nn = 10 %, exhibiting the 
increasing behavior with the channel energy. This result 
indicates that the disorder effect on conductance becomes 
more significant for low energies [Figs. 3(c) and 3(d)]. 
In addition, for a given channel energy, the localization 
length tends to decrease with increasing Uh- Thus, it 
is expected that a metal-to-insulator transition occurs 
as the adsorbate density increases. The charge density 
plots for the conducting channels also show the charac- 
teristics of localization. For hydrogenated graphene with 
L = 14 nm, the charge densities of the states around 
E = —0.6 eV are compared for different adsorbate den- 
sities in Figs. 4(a) and 4(b). For the low density of 1 
%, the energy states near E = —0.6 eV exhibits the ex- 
tended charge densities over the whole sample. Thus, 
this channel has the metallic conduction, with the local- 
ization length much larger than the sample size. As nn 
increases to 10 %, the localization is extended to higher 
energy states, reducing the localization length to a few 
nanometers. The localized behavior of the energy states 
around E — —0.6 eV is clearly seen in the plot of charge 
densities. 

From the results for (lng) in Fig. 3, the scaling func- 
tion /3(g) is drawn as a function of (lng) in Fig. 4(c). The 
calculated values for /3(g) for different adsorbate concen- 
trations and different energies are well fitted by an one- 
parameter function. For (lng) > 1, /3(g) is nearly zero, 
implying that graphene is in the 2D metallic phase, in 
which conductance is almost independent of the sample 
size. For (lng) < —2, as /3(g) is linearly proportional 
to (lng) with a slope of 1, the conductance follows the 
localization function, g ~ exp(— L/£), in the insulating 
phase. In the intermediate region, —2 < (lng) < 1, /3(g) 
is smoothly connected from the metallic to insulating 
phase, satisfying the hypothesis of the scaling function. 

We point out that /3(g) is always negative in Fig. 4(c), 
while its value is very small for (lng) > 0. Thus, the con- 
ductance tends to decrease continuously with increasing 
L, suggesting that an infinitely large graphene with H 
adorbates would act as an insulator. Note that our scal- 
ing function for short-range disorders is very different 
from those obtained in graphene subject to long-range 
disorders while it is similar to that derived for strong 
long-range impurities. 9 In the presence of long-range po- 
tential, the intervalley scattering of the Dirac fermions is 
generally suppressed, resulting in the positive beta func- 
tion which increases with decreasing g whereas converges 
to zero for large g. Thus, the conductance increases with 
sample size, robust against disorder, because none of the 
states can be localized. In contrast, hydrogen adsorbates 
in graphene act as short-range scatters and manifest the 
intervalley scattering. Due to significant backscattering, 
the Dirac fermion states are localized, giving vanishing 
conductance for a large system. 

In summary, our study shows that, even with low hy- 
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FIG. 4: (Color online) The charge densities (in clouds) of 
the states around E = —0.6 eV from the neutrality point are 
drawn for samples with the H concentrations of (a) 1 and (b) 
10 % and L — 14 nm. Dots stand for the positions of random 
H atoms, (c) The beta function is plotted as a function of 
(hxg) for different energies and different H concentrations. 



drogen concentrations of 5 — 10 %, a metal-insulator 
transition can occur by the localization of electron 
states. We find that the conductance decays expo- 
nentially with increasing of the sample size, satisfying 
the one-parameter scaling function in the 2D localiza- 
tion theory. This result provides an explanation for 
the recent observation of a metal-insulator transition on 
graphene-terminated SiC(OOlO) surface decorated with 
small amounts of atomic hydrogen/^ If the adsorbate 
density increases higher, the gap opening is likely to 
be a major cause for the insulating behavior, as ob- 
served by experiments! 12 ' 15 There is also the possibil- 
ity that adsorbate atoms are clustered at high densities, 
modifying the Dirac fermion nature into an insulating 
phase. Finally we point out that other adsorbates such 
as F atoms and CH3, C2H5, CH2OH, and OH molecules 
may also induce a metal-insulator transition, which sat- 
isfies the beta function in Fig. 4(c). Recent theoretical 
calculations 2 ^ showed that CH 3 , C 2 H 5 , and CH 2 OH 
molecules on graphene give rise to localized states in the 
midgap, similar to that of an H adsorbate, and low con- 
ductances around the Dirac point. Similarly, adsorbates 
such as F and OH directly interact with the host atoms, 
forming localized levels. As these defect levels are rather 
dispersive, there may be an asymmetry in conductance 
with respect to the neutrality point. 
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